Scanning Superfluid- Turbulence Cascade by Its Low- Temperature Cutoff 
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On the basis of recently proposed scenario of the transformation of the Kolmogorov cascade into 
the Kelvin-wave cascade, we develop a theory of low-temperature cutoff. The theory predicts a 
specific behavior of the quantized vortex line density, L, controlled by the frictional coefficient, 
ct{T) <^ 1, responsible for the cutoff. The curve \nL{\na) is found to directly reflect the structure 
of the cascade, revealing four qualitatively distinct wavenumber regions. Excellent agreement with 
recent experiment by Walmsley et al. [Phys. Rev. Lett. 99, 265302 (2007)] — in which L{T) has been 
measured down to T ~ 0.08 K — implies that the scenario of low-temperature superfluid turbulence 
is now experimentally validated, and allows to quantify the Kelvin-wave cascade spectrum. 
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In the studies of superfluid turbulence (ST) (for intro- 
duction, see, e.g., P, Q]), the case of very low tempera- 
tures remains intriguing and challenging Q . A wealth of 
theoretical work has been devoted to ST in this regime 
in thepast decade and especially in the past few years 
i, i, H 0, i, i, [M HH, [ll E [ll , and a detailed theo- 
retical picture seems to emerge. In this picture, the most 
interesting physics is associated with the short length 
scales, where quantized nature of vortices manifests it- 
self and the key role is played by Kelvin waves (KWs) — 
the distortion waves on quantized vortex lines. However, 
the feasibility of experimental confirmation of these pre- 
dictions with the current technique seemed vague, which 
to many made the theory appear purely academic. Until 
very recently, the only experimentally confirmed fact had 
been the existence of cascades in the T ^ limit [l^, UB] ■ 

In their recent remarkable work [l3|, Walmsley et al. 
measured the vortex line density L (the length of vor- 
tex lines per unit volume), of the Kolmogorov cascade in 
superfluid "^He down to temperatures T ^ 0.08 K. They 
found a striking temperature dependence of L and ar- 
gued that it should be due to the Kelvin-wave cascade 
E i, i, 0, B i, M MMMM being extended to 
higher wavenumbers with decreasing temperature. Bet- 
ter understanding of this process was formulated as a 
challenge for theorists. 

In this Letter, we meet the above-mentioned challenge. 
Recently, the two of us proposed a detailed scenario 
of how the zero-temperature Kolmogorov cascade trans- 
forms into the pure KW cascade. We argued that the 
transformation involves a chain of three intermediate 
regimes and thus occupies a significant interval in the 
wavenumber space. On the basis of this zero-temperature 
scenario, we now develop a theory of the low-temperature 
dissipative cutoff of the cascade. We show that the tem- 
perature dependence of L comes from the dependence of 
the wavelength scale Acuttofi , at which the cascade ceases 
due to the mutual friction of vortex lines with the nor- 
mal component, on the dimensionless friction coefficient 
a (xT^ (r — > 0). A very characteristic shape of the func- 
tion In L (In a) following from our analysis is in an excel- 



lent agreement with the experimental results of Walmsley 
et al. Remarkably, the shape of the curve In L(ln a) di- 
rectly reflects four qualitatively distinct wavenumber re- 
gions of the cascade. The agreement between our theory 
and the data of Ref. [l3] means that the scenario of low- 
temperature superfluid turbulence is now experimentally 
validated. Moreover, fitting the experimental data with 
theoretical curves allows us to quantify the KW cascade 
spectrum. 

A cascade of energy — the decay mechanism of both 
classical and superfluid turbulence — implies the existence 
of an inertial range — a significant range of length scales, 
Aon > A > Acutoff, in which all the dissipative mecha- 
nisms are weak rendering the system essentially conser- 
vative. The long- and the short-wavelength ends of the 
inertial range play a special role: almost all the energy of 
the system is concentrated at Aon, while the dissipation 
becomes appreciable only at A < Acutoff = Aoutoff(?^)- 
The relaxation process must lead to a transfer of the tur- 
bulent energy towards AoutofT, where it can be dissipated 
into heat. For the cascade regime to set in, it is also nec- 
essary that the dynamics within the inertial range is such 
that the "coUisional" kinetic time Toon(A), i.e. the time 
between elementary events of energy exchange at a cer- 
tain scale A, gets progressively shorter down the scales. 
In this case the decay is governed by the slowest kinetics 
at Aen, where the energy flux (per unit mass) e is formed, 
while the faster kinetic processes at shorter scales are able 
to instantly adjust to this flux supporting the transfer of 
energy between neighboring scales towards Acutoff - Thus, 
the cascade is a (quasi-)stationary regime in which the 
energy flux e is constant through the length scales, while 
the variation of e in time happens on the longest time 
scale of the system rooii(Aon)- 



In Ref. we argued that at the effective zero tem- 
perature the inertial range is highly nontrivial, since the 
energy transport mechanisms change several times within 
it. We first briefly outline the T = case (see IJ] for 
details), and then proceed to include the effect of mu- 
tual friction due to the flnite T. In this scenario, the key 
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dynamical parameter is 



A 



ln(/o/ao) > 1 , 



where oq is the vortex core radius and Iq is the typical 
separation between the vortex lines determined by the 
Kolmogorov energy flux e |14| . The parameter A controls 
the competition between the local self-induced motion 
of the vortex lines and the coupling between them. At 
the largest length scales in the inertial range, A ^ rg ~ 
A^^'^Iq, the vortex lines are coupled in bundles that move 
coherently mimicking the velocity profile of the classical- 
fluid Kolmogorov turbulence. (The indistinguishability 
of ST from its classical counterpart at large scales has 
been attracting much attention since it was first pre- 
dicted in Ref. [1].) At the scale tq, the self-induced 
motion of vortex lines becomes appreciable and the cas- 
cade enters the quantized regime, where each vortex line 
evolves independently (apart from reconnections). Be- 
fore the purely non-linear KW cascade [l^ can develop 
on separate vortex hues at fc ^ /c* ~ A^/^/?o the cascade 
experiences a complex transformation stage necessary to 
crossover from spatially organized bundles at k <^ r^^ 
and to generate the KWs. In the crossover, which (given 
A ~ 10) extends for about a decade in the wavenum- 
ber space, the vortex-line reconnections Q become the 
main driving mechanism of the energy transport in the 
cascade. Remarkably, the crossover range is further split 
into three (sub-)regimes distinguished by their specific 
types of reconnections: (1) reconnections of vortex-line 
bundles (r^ ^ <C <C fcb ~ l/A^/^/p); (2) reconnections 
between nearest-neighbor lines (fcb k <^ kc ^ A^/^/li^)] 
(3) self-reconnections on single lines {k^ k ^ k^,). It 
is important here that each regime features a distinct 
spectrum of KW amplitude bk summarized in Fig. [T] 

If the vortex lines were smooth, the vortex-line density 
L would be simply related to the interline separation as 
L = l^"^ . The increase of L due to the creation of KWs 
on the vortex lines is related to their spectrum by Q 



ln[L(a)/Lo] - 



ott(a) 



{hkkf dk/k 



(2) 



Here fccutoff ~ 1/Acutoff, k is the smallest wavenumber of 
the KW cascade (not to be confused with the smallest 
wavenumber of the Kolmogorov cascade) at which the 
concept of a definite cutoff scale is meaningful, and Lq is 
the "background" line density corresponding to kcutoS ^ 
k. There is an ambiguity in the definition of bk associated 
with the choice of the spectral width of the scale k, which 
is fixed in Eq. Q by setting the proportionality constant 
between the l.h.s. and the r.h.s. to unity. 

At r = 0, the cascade is cut off by the radiation of 
sound (at least in ''He) @ at the length scale Acutoff = Aph 
[ll[, As we show below, changing the temperature 
one controls Acutoff (7^) > Aph in Eq. ([2]), which allows one 
to scan the KW cascade observing qualitative changes in 
L(T) as Acutoff traverses different cascade regimes. The 
existence of a well-defined cutoff is due to the fact that 



the cascade is supported by rare kinetic events in the 
sense that Tcoh = Tcoii(e, k) is much larger than the KW 
oscillation period, Tpcr = Tpor(^)- The dissipative time 

'^dis = '^dis 

(a, k) ^ Tpor/a, as we show below, is the typ- 
ical time of the frictional decay of a KW at the scale k. 
Thus, the cutoff condition is 



Tdis(a, k) 



ou{£, k) , 



(3) 



which implies that the energy dissipation rate at a given 
wavenumber scale becomes comparable to the energy be- 
ing transferred to higher wavenumber scales per unit time 
by the cascade. It is this condition that defines the cut- 
off wavenumber fccutoff = ^cutoff (e, a). Decreasing T and 
thus a(r), one gradually increases fccutoff thereby scan- 
ning the cascade. In view of Eq. this in principle 
allows one to extract the KW spectrum. 




FIG. 1: Spectrum of Kelvin waves in the quantized regime 
as predicted in Ref. [ij and quantified here, apart from the 
regime (1), by the fit to experimental data, Fig. (2] In view of 
Eq. ((2)1, the constants C and ^2 can be found only in the com- 
bination « 0.049. The inertial range consists of a chain 
of cascades driven by different mechanisms: (1) reconnections 
of vortex-line bundles, (2) reconnections between nearest- 
neighbor vortex lines in a bundle, (3) self-reconnections on 
single vortex lines, (4) non-linear dynamics of single vortex 
lines without reconnections. The regimes (3) and (4) are fa- 
miliar in the context of non-structured vortex tangle decay 

aia. 

At finite T, dissipative dynamics of a vortex line el- 
ement is described by the equation (omitting the third 
term in the r.h.s., which is irrelevant for dissipation) [ll,0| 



s = v(s) + as' X [vn(s) — v(s)] 



(4) 



Here v(r) is the superfluid velocity field, Vn(r) is the nor- 
mal velocity field, s = s(^, t) is the time-evolving radius- 
vector of the vortex line element parameterized by the 
arc length, the dot and the prime denote differentiation 
with respect to time and the arc length, respectively. 

At a 1 the superfluid and normal components are 
strongly coupled and the KWs are suppressed. In this 
case, the cascade must cease before it enters the quantized 
regime, i.e. Acuttoff ^ ''o- In this Letter, we do not 
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discuss the strongly coupled case, since it has no control 
parameter for a theory. Consideration regarding energy 
spectra in this regime can be found in Ref. ■ 

If, however, the mutual friction is small, then S> 1 
gives the characteristic number of KW oscillations re- 
quired for the wave to decay. Indeed, to the first ap- 
proximation in 1/A, v(s) in Eq. ^ is given by the local 
induction approximation (LIA) [1], 



v(s) 



A;^ — s' X s" 
47r 



Ar = ln(i?/ao) 



(5) 



where k is the circulation quantum and R is the typical 
curvature radius. Taking into account that A^; is a very 
weak function of R, we shall treat it as a constant of the 
typical value A^ ~ A. As long as a <C 1, the disturbance 
of the normal component caused by vortex line motion 
can be neglected in Eq. (jlj , so that Vn (s) is irrelevant for 
KW dissipation. In this case, Eqs. (|3]), ([5]) give the rate 
at which the amplitude of a KW with wavenumber k 
decays due to the mutual friction: 



aujkbk 



(6) 



while the KW dispersion is ujk = (K/47r)Afc^. Here and 
below we omit factors of order unity, which are subject 
to the definition of the spectral width of the wave. These 
factors can not be found within our theory, but can be 
extracted from experimental data (as we do it below) and 
from numerical simulations. Since the energy per unit 
line length associated with the wave is Ek ^ Kpujkbf,, the 
power dissipated (per unit line length) at the scale ^ k^^ 
is given by 



(7) 



where p is the fluid density. In the following, we analyze 
ST as a{T) scans through the regimes shown in Fig.[TJ 

Regime (1). As it was already mentioned, at this stage 
the vortex lines are organized in bundles; the amplitudes 
of waves on these lines are given by bk ^ r^^k^^. 
One can estimate the total power lost due to the friction 
of these bundles at the wavenumber scale r^^ <^ k k\y 
per unit mass of the fluid as 



edis(A:) 



Pbl 



bj 



n(fc) 



a K 



A/ It- 



(8) 



Here, the first factor in the r.h.s. is associated with the 
correlation volume at this scale ~ b^k~^ [ij] and the 
second one stands for the number of vortex lines in the 
volume. Note that the dissipated power is constant at 
all the length scales within this regime and, since accord- 
ing to Ref. [31 the interline separation is related to the 
Kolmogorov flux hy e ^ Ak'^/Iq, it is simply given by 
ae. Thus, when a 1 the kinetic channel in the whole 
regime (1) becomes efficient and the cascade reaches the 
scale Ab , where the notion of bundles becomes meaning- 
less. That is the regime (1), as opposed to the regimes 
(2)- (4), is not actually scanned by a. Basically, the iner- 
tial range of the regime (1) develops as a whole while a 



evolves from the values of order unity to the values much 
smaller than unity. 

Regime (2). In the range k\, <^ k <^ kc , the spectrum 
of KWs is given by 6fc h{k\,/kf/^, which for the 
dissipated power yields 



edis{k) 



1 

'pit 



^ii{k) 



'A"^k^/lo 



(9) 



The condition edis(fc) ^ £ gives the cutoff wavenumber 

fccutoff = 6 A-^l^a-^'"^ (27r//o) , k^ < fceutoff < K , 

(10) 

where ^2 is some constant of order unity. Then, from 
Eq. ^ we obtain 



In 



i(fccutoff) 



C2(fcb/o)'[fccutofi7A:b-l] 



(11) 



Here, we set bk = C lo{kh/ky^^, where C is a constant 
of order unity. The overall magnitude of bk in the other 
regimes follows then from the continuity. 

Regime (3). In this regime, supported by self- 
reconnections, the spectrum is given by bk ^ k^^ (up 
to a logarithmic pre-factor). The corresponding energy 
balance condition yields the cutoff scale: 

fceutoff = 6 (Aa)"^/^ (27r/Zo), fcc < fceutoff < K ■ 

(12) 

With the logarithmic pre-factor taken into account 
[^], the spectrum in this regime reads bk = C[l + 
c§ ln(A;/fcc)]~"^/^(-\/fccfcb/fc)^o, where C3 is a constant of 
order unity. Then, the relative increase of vortex line 
density through this regime is given by 



-^(fceutoff) 
i(fcc) 



1 



In 



^'Cutoff 



(13) 



where v — C'^kck^ll/ c\. 

Regime (4). Since the spectrum of the purely nonlinear 
regime, bk ~ k~^^^, is steeper than the marginal bk ^ 
k^^ meaning that the integral in Eq. ^ builds up at 
the lower limit, as soon as fceutoff ^ fc* the line density 
^(fceutoff) starts to saturate and becomes independent of 
fceutoff at fceutoff ^ fc*- The energy balance gives the 
dependence fceutoff on a in the form 



fceutoff =e4A~3/4a-5/8(2^/Zo), k, 



cutoff 



> fc* 



(14) 



where ^4 is an unknown constant. The coeffi- 
cient in the KW spectrum is fixed by continuity 
with the previous regime, yielding bk — C[l + 

ciln(/c,//cc)]"^/^Vttfc*^^fc~^/^?o- Eq. 1^ thus yields 



cutoff ) 



(5C^ /2)fcefcb^^ 

1 + C3 ln(fc*/fcc) 



"-cutoff 



2/5' 



(15) 

The continuity of fceutoff leads to the following con- 
straints on the free coefficients in Eqs. PH)) . (I14p . 
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6 = J2 A'/^aJ'^, U = 6 AV^cyV^s, where 
ttc ~ A~^/^ is the value of the friction coefRcient at which 
the cascade is cut off at the crossover between the regimes 
(2) and (3) and a* ~ 1/A^ corresponds to the one be- 
tween (3) and (4). Introducing ctf, < 1, corresponding 
to the crossover from the regime (1) to (2), we rewrite 
Eqs. dm), dni), (HSl) in terms of a: 




a 



FIG. 2; (Color online.) Fit of the experimental data (squares 
and diamonds) adapted from Ref. {lo ~ 5x 10~^ cm giving 
A « 13). The high temperature measurements of Refs. [SJl, 
[2^ are represented by triangles and inverted triangles re- 
spectively. The form of a{T) is taken according to Ref. [31 
at r > 0.5K (roton scattering) and according to Ref [23] at 
T < 0.5K (phonon scattering). The fitting parameters are 
A2 « 0.36, A4, « 0.33, C3 ~ 2.9, u 0.18 with at « 2 x 10'^, 
ac ~ 2.5 X 10"*, Q, ~ 2 X 10"^ 

In conclusion, we discuss the qualitative form of 



L{a) in connection with the alternative scenario of the 
crossover from the Kolmogorov to the KW cascade pro- 
posed by L'vov et al. [l^. In Ref. 14!] we argued that this 
scenario mistakenly leaves out the reconnection-driven 
regimes. In this scenario the classical regime extends 
down to the scale of interline separation Iq, where it is 
immediately followed by the purely non-linear KW cas- 
cade [regime (4) in our notation]. Since the kinetics of 
the regime (4) are too weak to support the flux e at this 
scale, the vorticity was proposed in [13] to build up in 
the classical regime at the scales adjacent to Iq — the so- 
called bottleneck effect. Although this picture would also 
naturally predict a significant increase of L, the peculiar 
behavior of L{a) observed in [13] essentially rules out 
the scenario of Ref. [13] . The reason is that the dramatic 
rise of L happens only at a as low as 10"'^, while the bot- 
tleneck accumulation of classical vorticity must manifest 
itself as soon as KWs become an unavoidable relaxational 
channel, i.e. already at a < 1. 



Instead, L{a) exhibits the qualitative form peculiar to 
our scenario of Ref. [l3]. As a decreases from a 1 to 
the values significantly smaller than unity, the line den- 
sity L increases only by some factor close to unity (~ 1.5 
in the experiment), which reflects the formation of the 
regime (1) driven by the reconnections of vortex bun- 
dles. During the crossover from the region (1) to (2), the 
increase of L is minimal — a shoulder in the curve L{a) 
arises, in contrast to a large increase of L expected in the 
bottleneck picture. It is only well inside the region (2) 
that the increase of L becomes progressively pronounced, 
and, at the crossover to the region (3), the function L{a) 
achieves its maximal slope, determined by the fractaliza- 
tion of the vortex lines necessary to support the cascade 
within the interval (3). As the cutoff moves along the 
interval (3) towards higher wavenumbers, the slope of 
L{a) becomes less steep due to the decrease of the char- 
acteristic amplitude of KW turbulence. When the cutoff 
passes the crossover to the regime (4), the curve L{a) 
gradually levels. Fitting the experimental data fixes the 
values of all the dimensionless parameters, thus revealing 
the quantitative form of the KW spectrum. Fig. [T] Note, 
however, that in view of the fact that the experimental 
A ~ 10 is not so large we would expect certain systematic 
deviations, but a significant scatter of the experimental 
data does not allow us to assess them. 
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